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Abstract. We study the quantum energy of the Z-string in 2+1 dimensions 
using the phase shift formalism. Our main interest is the question of stability of 
' a Z-string carrying a finite fermion number. 

(N 

^ , PACS numbers: 03.65Sq, 03.70+k, 05.45.Yv, 11.27.-^d 

1. Introduction and motivation 

o 

• Z-strings were first discovered as solutions of the classical field equations of the 

, electroweak model by Nambu ^ in the context of bound pairs of magnetic monopoles. 

' Later on they were rediscovered - as independent objects in their own right - by 

■ ■ ' Vachaspati 0. 

^ • The main point under investigation in our study of Z-strings is their stability. If they 

D I are stable, they would be relevant for a variety of reasons: first of all, they would be 

i-S^ ' the first solitonic objects in the Standard Model to be found; given the importance and 

ubiquitousness of solitonic objects in effective field theories in general it is surprising 
that they seem to play no role in the Standard Model. A second observation that 
^ , might make them relevant is an alternative scenario of (electroweak) baryogenesis 

' proposed by Brandenberger et al P). The presence of networks made from Z-strings 

would make the requirement (for baryogenesis to happen at the electroweak transition) 
of a first-order electroweak phase transition obsolete. This is an attractive scenario 
since the electroweak transition is known not to be of first order. A third reason for 
studying their stability is that - since Z-strings end in magnetic monopoles - they might 
contribute to the primordial magnetic field. For a general overview of applications and 
properties of Z-strings along with a large collection of references cf. |4| . 

The structure of this paper is as follows: in section [21 we briefiy discuss the 
notion of stability in the presence of a conserved quantum number. Then, in 
section 121 we discuss the model under consideration, our method for computing the 
fermion determinant and some thoughts necessary to choose parameters properly in 
the 'D—2+1 dimensional theory. In section 01 we present the gauge and Higgs field 
ansatze used for computing the fermion determinant. Then in section we present 
the (preliminary) results available at the time of the QFEXT'05 conference. In section 
Owe present an outlook and some (preliminary) conclusions. 
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2. Stability 

The question of stabihty is of prime importance in gauging the importance of Z-strings. 
At this point one should keep in mind that the Standard Model does not provide a 
topological stabilization mechanism for string-like objects, unlike e.g. the Abelian 
Higgs model. Hence every extended object has to be stable on energetic grounds. In 
the purely bosonic, classical sector of the electroweak model, Z-strings are solutions 
of the classical equations of motion. But these solutions represent only a saddle point 
of the classical energy functional and not a minimum. Hence, the Z-strings can decay, 
e.g. by condensation of or W bosons along the string - so long as the value of the 
weak angle is close enough to its physical value; in the unphysical region sin > 0.9 
the Z-string actually is - classically - stable. These issues are discussed in depth and 
detail in the excellent review 0] . If one considers fermions in additions to the bosonic 
sector of the electroweak model, one finds that the Z-string actually binds fermions 
to its core, some of them rather tightly @]. Hence, one can investigate a new kind of 
stability p] : one can compare the total energy of the Z-string plus N bound fermions 
to the energy of N free fermions. If it is less, one has certainly found an interesting 
object, since even if the configuration under investigation decays further, it cannot 
simply decay to the vacuum (as it might in the absence of occupied bound states) 
since the fermion number is conserved and it has already been established that a 
configuration exists with energy below N times the free fermion number. However, 
if one wants to consider the bound state energy of fermions one also has to take into 
account the fermion determinant, since it arises at the same order of an h (loop) 
expansion as the fermion bound state energy. Thus the quantity which we consider in 
this talk is the difference AE{N) between the Z-string energy including the energy of 
N bound fermions and the energy of N free fermions: 

N 

AE{N) = (Z-string) -I- i^^c (Z-string) + J^Hujf'- \-m). (1) 

i=i 

Here Edass denotes the classical (bosonic) energy, Eyac denotes the renormalized 
vacuum polarization energy originating from the fermion determinant including effects 
of the counterterms and uij'^' denote the bound state energies. The sum runs over the 
occupied levels only. The fermion number of the Z-string is N. 

3. Technical prerequisites 

The model we consider for our computations differs in some respects from the full 
electroweak model. First of all, we consider our fcrmionic weak iso-doublets to be 
degenerate in mass. This is the most serious of our simplifications and cannot easily 
be gotten rid of, since without the isospin symmetry the problem doesn't have enough 
symmetry to allow a partial wave decomposition which is at the heart of our approach 
to computing the fermion determinant! . As a second simplification, we drop the 
hypercharge field from our consideration and only consider the SU(2) gauge field. This 
seems to be an innocuous simplification and we have techniques to deal with the U(l) 
field. A third simplification that is used only to simplify the algebra is the restriction 
of our fermionic sector to one weak iso-doublet. The fourth and final simplification is 

I This channel decomposition also applies to the bound state part of the spectrum. In each channel 
we have a finite number of bound states which can be determined using e.g. a shooting method. 
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that we consider the model m D=2+l dimensions§. Obviously, calculations in D=2+l 
are much simpler than in D=3+l due to the simplified UV divergence structure. And 
since this is only an exploratory investigation - to see whether it worthwhile to do the 
vastly more complicated calculation in D=3+l - it seems sensible to keep things simple 
where possible. But if there was no connection between results in 'D—2+1 and D=3+l, 
this investigation could not fulfil its purpose. Fortunately there is some evidence that 
the behaviour of energies and energy densities of string-like objects in T)=2+l can be 
a good guide to the behaviour of energies (per unit length) and energy densities in 
D=3-|-l. In our investigation of electromagnetic flux tubes [SlIT] we have learned that 
- given that the same renormalization conditions are used in D=2-|-I and D=3-|-l - the 
renormalized quantum energies are indeed very similar in their functional dependence 
on widths and fluxes. We use as a working hypothesis that the same is true in the 
electroweak model. 

The techniques used for expressing the vacuum polarization energy (which is 
related to the fermion determinant by dividing out the time interval T for which 
the determinant is evaluated) in terms of phase shifts from an associated scattering 
problem have been described extensively in |H] and shall be outlined here only briefly. 
The vacuum polarization energy can be renormalized effectively by realizing that if 
one replaces the full phase shifts by their Born approximation (of n*'' order) one gets 
the same result as by restricting the full one-loop vacuum polarization energy to the 
sum of Feynman diagrams with n external insertions of the background fields. Hence, 
in D=2+l, the vacuum polarization energy is given by 



-l^(|^f|-m)-i [ dkiV^^^m)J2-^ 



N 



b.s. " M " n=l 

N 

+ Y.Efl + EcT, (2) 

n=l 

where w^-*' denotes the fermion boimd state energies|j, 5m the (full) phase shift in 

angular momentum channel M, 5^^^^ its n*'' Born approximation; E'^pj^ is the energy 
contribution computed from Feynman diagrams with n external legs and i?cT is the 
energy resulting from the counterterms ^. Note that both the k integral on the 
one hand and the sum of Feynman diagrams plus counterterms on the other hand 
are separately finite. This in particular distinguishes our investigation from the first 
computation of the fermion determinant in the background of a Z-string performed 
by Groves et al |10| where determinant and counterterms were individually divergent 
functions of the (proper-time) cut-off parameter and a finite result was only obtained 
by combining these two quantities. Alas, they are known only numerically, hence this 
procedure is numerically not stable. The other difference is that we focus on occupying 
bound states, whereas Groves et al were mainly concerned with computing the fermion 
determinant. 

A last topic that merits discussion here is the question on how to choose the 
parameters of our model. In 3-1-1 dimensions, the gauge, Yukawa, Higgs self coupling 
and vacuum expectation value of the Higgs field can straightforwardly expressed using 

§ Nonetheless we use Dirac four-spinors to describe the fermions in our theory. 

II The bound state energies are determined from the first-order form of the Dirac equation. 

^ The calculation of the vacuum polarization energy per unit length in D=3+l uses the same phase 

shifts, Born approximations and bound state energies but different kinematical factors Iq]. Of course 

also Feynman diagram and counterterm contributions are different. 
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the fermion mass, the Higgs mass, the tree level mass of the W boson (the W boson 
mass at one-loop level is a prediction based on the tree level mass because of our 
choice of renormalization conditions) and the Fermi coupling Gp. In 2+1 dimensions, 
the masses can be unambiguously re-used; however, it is not entirely obvious how to 
choose the Fermi coupling. It is ultimately related to a cross section - a concept that 
would need some translation into two spatial dimensions. 

The aim of our 2+1 dimensional calculation is to be a guide to a full 3+1 
dimensional calculation. For the vacuum polarization energy we have assured this 
by using the same renormalization conditions as we'd use in a D=3+l calculation. 

We now fix the parameters in such a way that this is also true for the classical 
energy: in D=2+l we compute an energy E'^^^^, but in D=3+l we compute an energy 
per unit length, E^j^l^/L. Hence, we require 

^3+1 

^cUils = fundamental length, (3) 

where the fundamental length is given by half the Compton wave length of the fermion 
we integrate out. This makes sense from a physical perspective, since the fermion 
integrated out sets the scale beyond which spatial structures cannot be resolved any 
more. Hence a volume of thickness half the Compton wave length of this fermion is 
'perceived' as a surface. Also in the QED flux tube computations a relative factor 
of 7r/m appeared between the D=2+l dimensional energies and D=3+l dimensional 
energies per unit length. This prescription now allows to express the model parameters 
in terms of the physical parameters of the theory in D = 3+1. The question on how 
to choose an appropriate Gi;- in D=2+l can thus be avoided. 



4. Z-strings 

We compute the vacuum polarization energy for Higgs and gauge field configurations 
of a specific form+: 

'P+ \ -«./ff(p)cosCi +/p(p) 

cl^o J \ /H(p)sinCie''^ 

= r2/G(p)sin2^i, 
P 

^W+ = fe-'^fG{p)sin2^i. (4) 



Here p denotes the two-dimensional radius, p = y^r^ + y^, and ip is the unit vector in 
azimuthal direction. The Higgs vacuum expectation value is given by v. The gauge 
coupling is given by g. The functions fH{p),fp{p) and /g(p) are profile functions 
and are discussed in the remainder of this section. The requirement of finite classical 
energy necessitates for p — > both — » and fa 0. For p ^ oo, /g ^ 1 and 
= + l/pP ~^ 1 are required. As is changed from | to the configuration 
is changed continuously from the Z-string to a purely scalar configuration without 
winding. Note that in this process the gauge invariant length of the Higgs field 
does not change since it is independent of ^i. The classical energy is a continuous 
function of which - for /p = - has a maximum at = n/2 and decreases 

+ In absence of the hypercharge field the Z-field reduces to the W^-field. VF+ is the conventional 
charged W field, (po denotes the neutral Higgs field and the charged Higgs field. 
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continuously as — > 0. This illustrates our earlier statement that the Z-string is 
classically unstable and can unwind without hitting a topological barrier. 

The ansatz shown in (0J is a subset of a more general ansatz called the sphaleron 
square, cf. also [TT] , 

For our numerical investigations we cannot deal with general profile functions, 
but rather need also ansatze for the functions that have the proper behaviour for 
small and large p: 

_^ -4- 

fH{p) = '^-e , fp{p) = ape , /g(/o) = 1 - e ™c3 . (5) 

Altogether we have four parameters (plus ^i): three widths whtWqtWp and one 
amplitude ap - the amplitudes for f^, fu are fixed by requirements of finite classical 
energy mentioned above. 

5. Results 

In this section we want to present a couple of preliminary results for classical, vacuum 
polarization and bound state energies. The plots have in common that we have fixed 
S,i — ■^j i.e. we present results for the Z-string configuration. 

Couplings and Higgs vacuum expectation value are determined by our choice of 
masses and the Fermi coupling. We choose for the fermion mass 170 GeV, for the 
Higgs mass 115 GeV, for the tree level W boson mass 80 GeV and for the Fermi 
coupling lO-^GeV^ 

When studying the plots the following point has to be kept in mind: since is 
fixed, we have a four-parameter numerical problem. The plots are only standard two- 
dimensional plots, hence energies can only be plotted as a function of a single variable. 
We have chosen NOT to fix the other parameter values but plot the energies of all the 
configurations that we have available. Hence for each observable there is not a curve 
but a band corresponding to the ranges of the variables not represented on the x axis 
of the plot. In figure^we show in the left panel the renormalized vacuum polarization 
energy (in units of the fermion mass) including Feynman diagrams and counterterm 
contributions as function of the width of the neutral Higgs field, wh- This seems to 
be the predominant dependence. The contents of the right panel are slightly more 
difficult to explain: in ^ only one part of AE{N) depends on the fermion number of 
the configuration, namely 

N 

AE'-^-{N)=J2H'-\-rn). (6) 

Since we look for a stable object we restrict our choice of bound states in © to the 
most strongly bound states. In the right panel of figure ^ we now plot AE^-''-{N) 
for N= 5, 10 and 20. For the lowest curve we occupy all available bound states, 
denoted in the following as AE^"- (alV). Thus points along this curve can and will 
correspond to different fermion numbers. Whereas the curves with fixed fermion 
number level off as wh increases this latter curve keeps decreasing as wh increases. 
This is due to an amazing proliferation of bound states. Whereas for wh ~ 2 there is 
only a handful of bound states, for wh ~ 6 there can be easily more than 80 bound 
states. The dependence of the number of bound states seems to be roughly quadratic 
- as Wh increases, both more and more angular momentum channels contain bound 
states and the number in each individual channel keeps increasing, too. More recent 
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Figure 1. The left panel shows the (completely renormalized) vacuum polarization 
energy as a function of the width wjj, the right panel shows the energy gained by filling 
bound state levels relative to the same number of free fermions, X/j— i^l"^?^ I ~ 
for N=5,10,20 and all available bound states. Note that this latter curve corresponds 
to comparing configurations with different fermion number. 



investigations, performed after QFEXT'05, show that this behaviour continues at least 
up to widths wh ~ 12. In tliis parameter regime wc have found configurations with 
around 500 bound states in various angular momentum channels. 

In figure |2] we plot the classical energy, the vacuum polarization energy and 
Ai?^-*- (all). This figure shows clearly the different orders of magnitude that are 
involved. Furthermore it can be seen clearly that the effect of populating bound 
states by far outweighs the increase in energy due to taking into account the fermion 
determinant. 

When combining the classical energy with the vacuum polarization energy and 
AE''-''-{N) to form AE{N), one has to keep in mind that the fermions - in contrast 
to the bosons in this model - carry a colour quantum number and that fermions with 
different colour are energetically degenerate. Hence, both E^ac and AE^-''-{N) have 
to be multiplied by the number of colours Nc before they are added to Edass- The 
fermion number under consideration then is actually N x Nc- 

In figure 13 we use Nc = 9. In the left panel of figure El we have filled the 10 most 
strongly bound states. We find a minimum, but since we plot AE this minimum has 
to be below zero to indicate a stable object. So for both = 10 and = 30 (panel 
in the middle) the object under consideration is not stable. The situation is different 
for N = 50, since the minimum there is clearly below zero. Since N does not include 
the colour degeneracy this stable object actually carries fermion number 450. 

Also, a certain pattern seems to emerge from figure |21 as increases, the 
minimum of AE{N) moves towards larger values of wh and the value of AE{N) 
at the minimum decreases. This is of course due to the fact that the larger wh the 
more actual bound states are available and hence the possible gain in energy by filling 
these bound states also increases. It may even be possible that for sufficiently large wh 
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Figure 2. The figure shows the classical energies (upper band of points), vacuum 
polarization energies (denoted 'one-loop contributions', medium band of points) and 
Ai?'''^' (all) (lower band of points) as a function of the width wu- 
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Figure 3. In this figure we present our results for AE{N) as a function of the width 
wh for different values of N, from left to right N = 10,30,50. The actual fermion 
numbers are (again from left to right) 90,270,450. 



the energy gain is large enough to allow a stable object even for Nc = 3 to exist, but 
this is a question under current investigation and cannot be answered at the moment. 

A different investigation - results will have to be reported elsewhere - considers a 
heavy fermion with masses around l.bTeV instead of the top quark mass used here. 
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6. Conclusions 

For this contribution only preliminary data were available. Nevertheless we can state 
that - given a sufficient number of colours - we have found a very interesting object 
that is stable in the sense of section |21 It is very large and carries a gigantic fermion 
number (450). At the moment it is not clear what happens if we increase the size 
of the object further - will wc find so many bound states that maybe it is possible 
to find a stable object even for Nc = 3? Therefore a in-depth investigation of the 
parameter space is urgently needed, and is already under way |12j . Also, the great 
effort to investigate the D=3-|-l case is now fully justified. 
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